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ABSTRACT 

We study the formation of photospheric emission lines in O stars and show that the rectangular 
profiles, sometimes double peaked, that are observed for some stars are a direct consequence 
of rotation, and it is unnecessary to invoke an enhanced density structure in the equatorial re- 
gions. Emission lines, such as N IV A4058 and the N III AA4634-4640-4642 multiplet, exhibit 
non-standard "limb darkening" laws. The lines can be in absorption for rays striking the center 
of the star and in emission for rays near the limb. Weak features in the flux spectrum do not 
necessarily indicate an intrinsically weak feature - instead the feature can be weak because of 
cancellation between absorption in "core" rays and emission from rays near the limb. Rotation 
also modifies line profiles of wind diagnostics such as He II A4686 and Ha and should not be 
neglected when inferring the actual stratification, level and nature of wind structures. 
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1 INTRODUCTION 

Many O stars exhibit rapid rotation that broadens their line profiles. 
In the absence of a wind, the influence of rotation can be accurately 
computed by integrating the intensities from the stellar disk with al- 
lowance for the projected velocity in the direction of the observer. 
For convenience, and speed, it is customary to allow for the ef- 
fects of rotation by convolving the flux spectrum with a rotational 
broadening function that accounts for limb darkening. This proce- 
dure, which assumes that the line profile does not change across 
the disk (i.e., t he continuu m and line have the same limb darken- 
ing laws) (e.g., iGravl 19921) , generally works very well. As our own 
tests have shown, the simple convolution procedure is sufficiently 
accurate for the majority of photospheric absorption profiles in O 
stars, although for the most precise work the disk integration pro- 
cedure should be used. 

In general, O stars exhibit a wind which modifies the photo- 
spheric Ha profile (possibly driving Ha into emission) and gener- 
ates numerous UV P Cygni profiles. Because the wind is extended, 
and because the rotation rate declines with radius, the simple con- 
volution technique is not valid. In such a case it is necessary to 
take the 2D structure of the wind into account. In the simplest case 
one can assume that rotation simply affects the velocity structure 
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of the wind while maintaining an almost spherical structure, while 
for rapid rotation it is expected that the density structure will also 
be altered. How the 2D density structure is altered is very com- 
plex, and depends on the rotation rate relative to the critical rota- 
tion rate, the rotation rate relative to the terminal wind velocity, 
and the closeness of the star to the Eddington limit. Depending on 
their values, it is possible to get either density enh ancements in the 
equatorial or polar flows (Maeder & Meynel 2000). A prolate wind 
structure is expected unless the number of lines (or more correctly 
the flux [line] mean opacity) increases sufficiently rapidly towards 
the equator to overcom e the decreasing equatoria l flux arising from 
gravity darkening (e.g- lMaeder & Mevnetll200oT) . 

In many O stars intense N III A4634-4642 emissio n is seen 
and is used as a criterion to define the so-called "f ' class (Walborn 
ll97lMSotaetal.l201 1 ). When resolved, these lines can show a rect- 
angular (or trapezoidal-like) profiltQ, with the stronger components 
exhibiting a double peaked struc t ure. In the spectra of O super- 
giants presented by iBouret et alj < |2012|) . ( Pup, HD 16691, and 
HD 210839 show broad, non-gaussian profiles for the N III multi- 
plet, although HD 15570, HD 163758 and HD 192639 exhibit more 
gaussian-like profiles. The mechanisms driving the N III lines into 



1 For simplicity we refer to these profiles as non-Gaussian, while other 
profiles will be referred to as Gaussian, even if they do not exhibit extended 
wings. 
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emissi on have been discussed previously bv lRivero Gonzalez et al.l 
l|). The dominant mechanism depends on the precise stellar 
param eters - both continuum fluo rescence (e.g., the Swings mech- 
anism, iBruccato & Mihalas|[l97ll) and dielectronic recombination 
(previously d iscussed for these lines by Mihalas et alj[T972l) can 
be important. iRivero Gonzalez et al also indicate that the 

strength of the lines can be influenced by interactions between the 
N III and O III resonance lines. 

In this paper we discuss the profiles of photospheric emission 
lines in O stars, and show that rotation can produce rectangular-like 
profiles that often exhibit double-peaked profiles. The phenomena 
is not restricted to the N III A4634-4642 complex — in £ PupN IV 
A4058 and Si IV AA4089,4116 also exhibit non-gaussian profiles. 
We show that the observed line profiles are a direct consequence 
of the lines being photospheric and being in emission. The lines do 
not show "classic" limb darkening - rather the intensity distribu- 
tion across the star can be flat, and may even show limb brighten- 
ing. In addition, we further examine the influence of rotation on the 
Ha and He II 4686 line profiles. Rotation produces observable ef- 
fects on the line profiles, and these effects are seen. Because of the 
influence of rotation on the wind velocity law (which is now 2D) 
wind line profiles potentially depend on both v and sin i rather than 
v sin i. To model these emission lines on O stars, it is absolutely es- 
sential to properly account for the effects of rotation in O stars that 
exhibit rapid rotation. It is especially true for rapid rotators such as 
£ Pup which has usini ~ 210 kms -1 . 



2 PREVIOUS WORK 

The influence of rotation on O star profil e s has been the subject of 
several different studies. IPetrenz & Puis! i ll 996!) studied the influ- 
ence of rotation on Ha. Their work was mainly concerned with the 
influence of rotation on the wind profile. Rotation directly influ- 
ences the velocity structure (the velocity along a given sight line 
is not necessarily monotonic) and can potentially alter the den- 
sity structure of the wind throug h direct means (as in the wind - 
compressed disk (WCD) model o f iBiorkman & Cassinellil dl993l) ) 
or indirectly by introducing a latitudinal variation of temperature 
and surface gravity which then causes a latitudinal variation of the 
mass loss rate and the velocity law. 

In the work o f |Petrenz& Pulsl d 19961) emission in Ha was sig- 
nificantly modified by the enhancement of density in the equato- 
rial plane arising from the WCD ef fect. However, this ef fect was 
likely substantially overestimated — lOwocki et alj l ll998h showed 
that non-radial forces will inhibit both the formation of a wind- 
compressed disk and a wind-compressed region. As a consequence, 
the most important influence of rotation on the wind arises from ro- 
tation introducing a latitudinal dependence of temperature and sur- 
face gravity, but this effect remains relatively small unless Vmt > 
O.SKrit. Two potentially important effects for wind lines are the de- 
crease of V4, with r which reduces the influence of rotational broad- 
ening on wind lines, and the resonance-zone effect which arises 
because rotation intro duces a non-monotoc ity in the velocity law 
alone some sightlines dPetrenz & Pulsl 19961) . Given the current un- 
certainties in our understanding of ID winds, it is fair to say that 
the structure of 2D winds, and the azimuthal variation of density 
and mass-loss, is still very uncertain. 

Another study o f the influence of rotation on line profiles in 
O stars was made by iBusche & Hillieri d2005l) . They examined a 
range of profiles, showed that the absorption profiles were ade- 
quately modeled via the convolution technique, and confirmed the 



important influence of rotation on Ha and He II A4686 line profiles. 
While the analysis was general, and applied to the whole spectrum, 
they did not study the formation of photospheric emission lines. 



3 TECHNIQUE 

To study the influen c e of r otation we use the 2D code developed 
by IBusche & HillieJ ( 120051) . This code computes the formal solu- 
tion (and hence line profiles) for axisymmetric 2D winds. It uses as 
input ID opacities and emis sivities computed with t he ID formal 
solution code, CMF_FLUX dBusche & Hillierll2005h . It has been 
parallelized, and runs efficiently on workstations with multiple pro- 
cessors. 

In the present calculation we assume that rotation does not 
influence the density structure — this assumption is a reasonable 
first approximation since the rotation rates are less than 50% of the 
critical rotation rate, and much less than the wind terminal veloc- 
ities. We also neglect the latitudinal variation of temperature and 
surface gravity with latitude — these are likely not to have a large 
effect since the rotational velocity of most O stars is generally a 
factor of 2 (or more ) less than the crit i cal ro tational velocity. Fol- 
lowing the work of IBusche & Hillie ri J2005h we assume that the 
star rotates as a solid body below 20kms~ 1 , and that angular mo- 
mentum about the center of the star is "conserved" at higher ve- 
lo cities (that is, V4, oc R (V = 20 km s -1 )/V). Following the work 
of lOwocki et al. I dl998h we assume non-radial forces inhibit disk 
formation and hence set the polar velocity component to zero. The 
last two assumptions only influence the formation of wind lines - 
they do not affect photospheric lines. For simplicity, we neglect the 
influence of rotation on the non-LTE po pulations - such effects are 
likely weak, and have been discussed bv lPetrenz & Pulsl d2000h . 



4 OBSERVATIONS AND MODEL 

To illustrate the influence of rotation we use £ Puppis which is a fast 
rotator with v sin i ~ 210km s -1 . This sta r has b e en extensively 
analyz ed, most recently by iBouret et al .1 d2012h . iNaiarro et al .1 
d201lh have also extensively analyzed this star and obtained sim- 
ila r results. The mode l parameters we adopt are from the paper 
of IBouret et al.l d2012h . and are as follows: T eff = 40 000 K, 
loggc = 3.64 (centrifugally corrected value of log g), M = 
2.0 x lO^Meyr" 1 , logL/L Q = 5.91, = 0.05, v cl = 
100 km s" 1 , N(He)/N(H)= 0.16, Z(C)= 2.86 x 10~ 4 , Z(N) = 
1.05 x 10~ 2 andZ(0)= 1.30 x 10~ 3 (Z denotes mass fraction)Q 
/oo and w c i control the clumping in the wind and have little influ- 
ence on the photospheric profiles. For illustration purposes, we use 
the same observational data (which covers 3800 to 8800 A at a res- 
olution of 48,000 and a signal-to-noise in excess of 100 to 1 per 
pixel) as presented in that paper. 

£ Puppis is an ideal candidate for the study. Its parameters 
are well known (except for uncertaintie s induced by the the uncer - 
tainty in its distance; see discussion in lSchilbach & R6seill2008I) . 
and from variability studies it is most likel y observed edge on (i.e ., 
rotation axis perpendicular to line of sight; iHowarth et all dl995h ). 
It exhibits several emission lines such as N IV A4058, N III A4634- 
4642, and Si IV AA4089,4116 that show non-gaussian profiles - 
generally rectangular but possibly exhibiting enhanced emission 



2 Expressed as number fractions: N(C)/N(H)= 
1.2 X 1(T 3 and N(0)/N(H)= 1.4 X 10~ 4 . 



4.0 X 10" 5 ,N(N)/N(H)= 
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near the extremes of the emission profile. Further, C IV 5801, 5812 
exhibits complex profiles while both Ha and He II are in emission. 

Observed profiles for a selection of lines in £ Pup are shown 
in Fig, [7] Also shown are the profiles obtained from a simple con- 
volution of the ID model flux with the rotation profile, and the pro- 
file computed from a full integration of the intensities across the 
disk of the star (hereafter disk profiles). For simplicity we adopted 
a macroturbulent velocitjf] of 0km s^ dBouret et alj|2012l adopted 
90 km s _1 ) and a rotation velocity of 240 km s _1 (and sini = 1) to 
better match the width of the rec tangular emission p rofiles rather 
than the 210 km s -1 adopted by iBouret et alj d2012l) . As readily 
apparent, the disk profiles generally show much better agreement 
with the observed profiles than do the profiles obtained through 
convolution. In particular, N IV A4058 exhibits a double peaked 
profile, as observed. Similarly, the N III lines show rectangular pro- 
files. Lines affected by the wind also show improvements in pro- 
file shapes — He II 4686 exhibits evidence for "absorption" on the 
blue side which is clearly seen in the observation but is not present 
in the profiles created by convolution. Profiles computed assuming 
a higher rotation velocity of 480 km s -1 and sin i — 0.5 so that 
v sin i = 240 are essentially identical to those computed using the 
lower rotation velocity. In reality, since the wind and star properties 
would change with such a high rotation rate, the profiles would not 
be identical. 

It is not surprising that the convolution procedure does not 
work well for Ha and He II A4686. Both lines are significantly in- 
fluenced by the wind, and line transfer in the wind will be signifi- 
cantly affected by the presence of rotation. 

From Fig, Q] is is apparent that predicted disk profiles do not 
agree precisely with those observed. This likely arises because of 
errors in the adopted stellar parameters and abundances, errors in 
the atomic data, and missing physics. A full discussion of the uncer- 
tainties associated w ith the parameters of ( Pup has been given by 
IBouret et alj d2012l) . As the purpose of this paper is to understand 
the physics giving rise to the observed rectangular profiles, and to 
re-examine the influence of rotation, we do not attempt to readjust 
model parameters. We note, for example, that the absorption asso- 
ciated with Ha is sensitive to the radial variation of clumping with 
radius dHillier et alj2003l;|Puls et alj2006l : lNaiarro et al.l201ll) and 



porosit y effects l Sundqvist et alj201lh . Our own work on O super- 
giants IBouret et al j|2012h also shows that it is model sensitive (it 
varies significantly when the stellar parameters are varied over their 
error range). 

To understand the observed emission line profile, and why 
convolution does not work we consider N IV A4058. In Fig. [2] we 
show the intensity variation at line center (no rotation) and in the 
adjacent continuum as a function of impact parameter. As appar- 
ent, the intensity variation in the continuum and the line is very 
different — the continuum shows "classic" limb darkening, while 
the intensity in N IV A4058 is almost flat. Since the "limb darken- 
ing laws" for the continuum and line differ substantially, we cannot 
use the simple convolution procedure to obtain the profile when the 
star is rotating. Rather, the observed line profile must be obtained 
by integrating the intensities across the disk. 

From Fig.[2]it is also apparent that the wind makes relatively 
little contribution to N IV A4058. N IV A4058 is a photospheric 
emission line, and is in emission because of non-LTE effects. In 



3 The use of a (constant) macroturbulent velocity is, of course, another 
approximation which will modify line profiles and hence influence conclu- 
sions that can be drawn from them. 
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Figure 2. Illustration of the "limb darkening law" for the N IV A4058 line 
at line center (blue curve) and the adjacent continuum (A4056; red curve). 
The continuum shows a classic limb-darkening law. Conversely, because of 
non-LTE effects, the N IV A4058 line intensity is almost constant across the 
star and brightens near the limb. 



Fig. [3] we illustrate the origin of the N IV A4058 as a function of 
the Rosseland mean optical depth, and velocity. This diagram also 
confirms that N IV A4058 is photospheric — the bulk of the line 
originates inside lOkms -1 . 

To better understand why N IV A4058 exhibits a non-gaussian 
profile we can examine the center-to-limb variation of the intensi- 
ties as a function of velocity, and this is shown in Fig.|4]for a range 
of impact parameters. As readily apparent, the N IV A4058 profile 
shows substantial variation with impact parameter - for the core 
ray the line is in absorption while it is strongly in emission near 
the limb. Since the observed flux is a complicated weighting of the 
intensities, we illustrate in Fig. [5] the growth in the spectrum with 
impact parameter. When integrated over the inner core rays (for all 
azimuthal angles) we see an absorption line spectrum, but as we 
extend the integration zone to higher impact parameters we begin 
to see emission, and eventually we get the double peaked profile, 
similar to what is observed. 

Another line that exhibits a rectangular-shaped profile is N III 
A4634. While the line is in emission at all impact parameters 
(Fig.[6j, it is much more in emission at the limb. It is this enhanced 
emission from "large" impact parameters that changes the profile 
shape from Gaussian to rectangular (Fig. [7j- 

To complete our study we illustrate two other lines. The first 
is C IV A5801 which shows an absorption component together with 
extended wind emission. In many O stars this li ne is very difficult to 
model, being sensitive to both line blanketing dBouret et alj|2012l) 
and model parameters. This line, like N IV A4058, is in absorp- 
tion for core rays (Fig. [§}, but in emission for rays near the limb. 
However, in contrast to N IV A4058, the line remains in absorption 
longer (i.e., to higher impact parameters) before switching to emis- 
sion. The central part of the line in the full spectrum thus remains in 
absorption (Fig.[9jl. Careful study of Fig.[9]reveals extended emis- 
sion line wings arising from the wind. Given the transition from 
absorption to emission as we move to higher impact parameters, 
and since what we observed is a weighted sum of this variation, it 
is not surprising that the line is model sensitive. 

From Fig. Q] we see that the C IV absorption profiles are nar- 
rower than what one would predict based of the rotation rate of 
the stars. Although not so evident in the current model, disk pro- 
files for C IV A5801 are generally narrower than those obtained 
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Figure 1. Illustration of a sample of line profiles. The red curve shows the spectrum of £ Pup, the blue curve shows the spectrum computed using the 
formal integral of the 2D model with vsini =240 kms -1 , and the green curve show profiles obtained through convolution of the ID-flux spectrum with 
v sin i = 240 km s" 1 . From left to right, top to bottom we show N IV A4058, Si IV AA4089,41 16, N III AA4634-4640-4642, He II A4686, C IV AA5081, 5812 
and Ha. In all cases the profiles computed with the 2D model show better agreement than the profiles obtained via convolution with that observed. In particular 
notice the distinct non-gaussian profiles associated with N IV A4058, Si IV A 4089, 4116, and N III AA4634-4640-4642. On this plot wavelengths are vacuum, 
but we use air wavelengths (as is the usual practice) when identifying lines. 



via convolution. To get a better match for the present model would 
require stronger C IV absorption for core rays, and (possibly) en- 
hanced emission for the limb rays. Because the absorption is pref- 
erentially confined to the core rays, the full disk absorption profiles 
would be narrower than for pure absorption lines. 

The final line studied is He II A4686. This line exhibits strong 
emission, but in £ Pup (and some other stars) it shows a notch on the 
blue side. As shown in Fig. [T] this notch can be explained by cor- 
rectly accounting for the effects of rotation. Like N IV A4058 and 



C IV A5801 this line also transitions from absorption for low impact 
parameters to emission for high impact parameters (Fig.|10|>. How- 
ever there are important distinctions in its behavior from these two 
lines. First, we see extended absorption to the blue, arising from the 
stellar wind. Second, wind emission is very important and is what 
drives the line into emission (Fig. II lb . 
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Figure 3. Illustration of the origin of the N IV A4058 line at line center (blue curve) and the adjacent continuum (A4057; red curve) in r and (radial) velocity 
space for a ID non-rotating star. Each data value (bin) provides the contribution to the total flux from that pixel - the total flux is the sum over all pixels. The 
two vertical lines in each figure indicate the location of tr oss = 2/3 (left line) and the location where the wind velocity matches the isothermal sound speed 
(~ 23kms~ 1 ). 
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Figure 4. Intensity along the equator as a function of impact parame- 
ter for N IV A4058. The rotation velocity is 240 km s . From top to 
bottom we show the intensity forp/R* = 0.0, 0.039, 0.733, 0.965, 
0.983, 0.997, 1.007, 1.0127, and 1.062 where R* is the radius of the 
star at a Rosseland mean optical depth of 2/3. Notice how the line is 
in absorption near the center of the star but goes into emission near 
the limb. We only show the curves for half the star - curves for the 
equatorial region rotating towards the observer are similar (although 
not identical) to those shown above. 
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Figure 5. Plot showing the spectrum computed by integrating from 
p = to p(i) where p(i)/R* = 0.039, 0.733, 0.965, 0.983, 
0.997, 1.007, 1.0127, 1.062 and 100 (=R ma x/Rstar). The curves for 
p(i)/R* = 0.039 (bottom) and p(i)/R* = 0.733 have been 
scaled upwards by factors of 7.5 and 1.7, respectively. These plots, 
together with Fig. [4] illustrate how the observed profile is a compli- 
cated blend of absorption and limb emission. 
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Figure 6. As for Fig. \4\ but illustrating N III A4634. The strength 
of the emission relative to the continuum is significantly enhanced 
towards the limb of the star. 
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Figure 7. As for Fig.|5]but illustrating N III A4634. At ~400 kms" 1 
we see a blend of two other N III lines (AA4641, 4642) that belong 
to the same multiplet. 
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Figure 8. As for Fig. [4] but illustrating C IV A5801. The strength 
of the emission relative to the continuum is significantly enhanced 
towards the limb of the star. 
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Figure 9. As for Fig. [5] but illustrating C IV A5801. At 
~500kms — 1 we see the other component (C IV A5812) of the dou- 
blet. In the top spectrum an additional broad emission component 
arising from the wind is evident. 
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Figure 10. As for Fig.|4]but illustrating He II A4686. Because of the 
wind, these profiles exhibit absorption in the blue. 




-1000 



-500 



500 



1000 





V(km s" 1 ) 

Figure 11. As for Fig.|5]but illustrating He II A4686. The additional 
curve (second from top) is for p/R t = 2.513, and was included 
because this line is significantly influenced by wind emission. 



5 CONCLUSION 

We have demonstrated that rotation alone can explain the rectangu- 
lar profiles, sometimes double-peaked, that are observed for some 
optical emission lines in O stars. It is not necessary to invoke depar- 
tures of the density structure from spherical symmetry to explain 
these profiles. To properly compute the emission line profiles, it 
is necessary to integrate the intensities across the rotating stellar 
disk - the classic convolution of the flux does not work because 
the line profiles vary substantially across the stellar disc (i.e., with 
impact parameter). In some cases the line is in absorption for rays 
striking the center of the star, while it is strongly in emission for 
rays near the limb. For N IV A4058 in £ Pup, the observed weak 
double peaked profile is actually a result of near perfect cancella- 
tion between absorption at low impact parameters, and emission at 
high impact parameters. This illustrates that the variation in inten- 
sity with impact parameter can provide additional insights into line 
formation. 

From thei r stud y of the N III lines in HD 16691, 
|Pe Becker etaLI ( 120091) concluded that the lines arise from close to 
the star in a large scale corotating structure. While we confirm that 
rotation is crucial in explaining the observed line shapes, we find 
that the N III lines are primely of photospheric origin in these su- 
pergiants and that rotation alone can explain their observed shapes. 



For most analyses it is adequate to use the normal convo- 
lution technique when modeling absorption lines. However, there 
are subtle differences between disk profiles, and the profiles ob- 
tained with convolution, and these could be important when us- 
ing absorption lines to study line asymmetries and the possible 
influence of the velocity field at the base of the wind. Following 
a suggestion by Francisco Najarro, and due to the intense inter- 
est with the mass of the most m assive stars, we also r an a model 
for R136a3, which was found bv lCrowther et alj bold) to have an 
initial mass of ~ 165 Mq and a v sin i ~ 200 km s" 1 . Profiles 
computed using a 2D model to take into account rotation showed 
only small differences from profiles computed using the convolu- 
tion tech nique, and are too smal l to significantly affect conclusions 
drawn bv lCrowther et alj fcOld) . 

Although already well known, we also highlighted the influ- 
ence rotation has on main optical wind diagnostics - Ha and He II 
A4686. While the correct allowance for rotation does not funda- 
mentally alter the line strength (and hence, for example, inferred 
mass-loss rates) it does significantly alter the line profile. This is of 
crucial importance - it means that Ha and He II A4686 emission 
line profiles in rapidly rotating O stars cannot be used to infer vari- 
ations in clumping with distance and the importance of velocity- 
porosity unless rotation is accurately taken into account. Further, 
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rotation will influence the accuracy of parameters, such as ft (in the 
classic wind-law v(r) = (1 — r/R^i) 13 ). 

The correct treatment of rotation is also important for UV res- 
onance lines formed in the wind. As photospheric lines can influ- 
ence the resonance line profiles, it is important to take rotation into 
account. However, a simple convolution overestimates the effect of 
rotation on the wind lines and can adversely affect estimates of the 
terminal velocity, and microturbulence within the wind. 

Since it is relatively easy to treat the gross influence of rotation 
accurately, this should become the norm in future studies of O stars. 
Issues related to the influence of rotation on the azimuthal density 
structure, and the dependence of clumping parameters on rotation, 
will require further studies. 
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